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Abstract 

The 'extended Aharonov-Bohm (AB) period' recently proposed by Kusak- 
abe and Aoki [J. Phys. Soc. Jpn 65, 2772 (1996)] is extensively studied 
numerically for finite size systems of strongly correlated electrons. While the 
extended AB period is the system length times the flux quantum for nonin- 
teracting systems, we have found the existence of the boundary across which 
the period is halved or another boundary into an even shorter period on the 
phase diagram for these models. If we compare this result with the phase di- 
agram predicted from the Tomonaga-Luttinger theory, devised for low-energy 
physics, the halved period (or shorter periods) has a one-to-one correspon- 
dence to the existence of the pairing (phase separation or metal-insulator 
transition) in these models. We have also found for the t-J model that 
the extended AB period does not change across the integrable-nonintegrable 
boundary despite the totally different level statistics. 
PACS numbers: 74.20.Mn, 71.10.Hf 
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I. INTRODUCTION 



Aharonov-Bohm (AB) effect is an interesting probe for various phenomena in condensed- 
matter physics, since this effect explores 'global' responses (over large changes in the gauge 
flux) rather than linear responses against an infinitesimal change. In this light, it should be 
especially interesting to look at strongly correlated systems such as the Hubbard model or the 
t-J model from the AB effect. A natural interest is how the occurrence of superconductivity, 
if any, should be reflected in the AB effect. The response of the system to the AB flux 
has long been recognized as the 'anomalous flux quantization' as a tool for detecting the 
Cooper pairing ever since Byers and Yang proposed the idea.S Let us pierce an AB flux, 
$, through the torus as shown in Figure.|l]. In a normal state, the ground-state energy is a 
periodic function of $ with period $o = hc/e with minima at $/<3>o = 0, ±1, ±2, ■ • •. In a 
superconducting phase, new stable states emerge at $/$o = ±2' -*-| ' ' ' ^ ue ^° ^ ne P a i rrn g; 
which makes the period $o/2. The method has been applied to several mo delfJHi. 

The response to an infinitesimal AB flux can also detect the Drude weight iM or the 
superfluid weight The Drude weight is a measure of the ratio of the density of mobile 

charge carriers to their mass, while the superfluid weight measures the ratio of the superfluid 
density to mass. For a clean system, we expect D = D s = for an insulator, D is finite and 
D s = for a metal, and D and D s are both finite for a superconductor. 

As a test for pairing, however, the conventional anomalous flux quantization can be 
ambiguous, since the method is, in some cases, not sensitive enough to determine phase 
boundaries, or because finite-size effects obscure the test. For example, Ferretti et al. shows 
that both repulsive and attractive Hubbard model behave as a superconductors in the sense 
of Byers- Yang arguments at T=0.i We can see this in the left corner of Fig.||]. 

In order to shed a new light, two of the present authors have recently proposed an 
entirely different method of detecting bound electron states (see Fig.0) .III The idea stemmed 
from Sutherland's study for one- dimensional (ID) Heisenberg magnets!! He looks at what 
happens to a many-body state when we twist the boundary condition as 
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0(.-.,x i + iV a ,...) = e i2 ^V(---,^v) ) 

for A^ a -site lattice shifts (or 'boosts') the total momentum of an iV particle state. For a 
noninteracting particles, the boost by $ = $o shifts each one-particle momentum by 27r /N a , 
so that $ = iV a $o will bring the set of k points back to the original position. 

For an interacting particles we can expect that $ = iV a $o will also recover the original 
eigenstate when all of the N particles move individually. By contrast, the response will be 
different if the particles form iV-bound states, since bound iV particles should have a well- 
defined center-of-mass momentum, whose response to the AB flux will be A" times greater. 
A twist of $ = N a &o/N will then suffice to bring the energy back to the initial value. 

For the Heisenberg spin chain Sutherland^! has in fact confirmed this behavior using 
the Bethe-ansatz analysis. Kusakabe and Aokii have then proposed that the idea for spin 
systems may be extended to electron systems, such as the Hubbard model, for detecting 
bound-electron states. In addition to the Bethe-ansatz analysis for electrons, the method 
is also numerically implemented in order to reveal the phase diagram of strongly correlated 
systems in (pseudo-)lD and possibly in 2D systems, where the Bethe ansatz is inapplicable 
but we can still keep track of the ^-dependence of the state over the 'extended zone' (0 < 
$ < A" a $ ) well beyond the one period (0 < $ < <3> )- 

Thus this technique is called 'the extended AB period method'. It is conceptually inter- 
esting to keep track of the states along the extended zone, since this questions an even more 
global response of a many-body system than the anomalous flux quantization. 

With this method a normal-to-superconducting transition is expected to be characterized 
as a halving of the extended period. For the ID Hubbard model, which is an integrable 
system, this has been shown indeed to be the case The spectral flows have also been 
analyzed in other integrable systems by Fukui and Kawakami in the context of the exclusion 
statistics, such as the Haldane-Shastry mo del0 and the supersymmetric t-J model with 
1/r 2 interaction.!^ On the other hand, the numerical implementation has been applied to 
non-integrable systems. @0 
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In the present paper, we investigate extensively the extended AB method by applying 
it to one- dimensional (1) t-J model, (2) t-J ladder model, (3) t-J-J' model, and (4) the 
extended Hubbard model to systematically look at the behavior of the extended period. 

We address ourselves two fundamental questions: 

(i) what is the difference between integrable and non-integrable systems in terms of the 
extended AB flow, 

(ii) can we actually detect normal-superconducting transition or phase separation in such 
nontrivial models as enumerated above, where a quantum phase transition (at T = 0) is 
believed to occur at a finite strength of the interaction. 

As for the first question, we can naively expect drastically different behaviors for the two 
classes of systems, since the extended AB period is dominated by whether the level crossings 
in the adiabatic flow line turn into anti-crossings, while the quantum adiabatic theorem dic- 
tates that the adiabatic evolution of a level is free from level crossings (except for accidental 
ones) when the system is nonintegrable. We can in fact characterize the avoided crossings 
in terms of the level statistics (distribution of the spacing of adjacent eigenenergies) , which 
is thought to be universally Gaussian orthogonal (unitary or symplectic depending on the 
symmetry of the Hamiltonian) for a nonintegrable system while an integrable system has 
Poisson's distribution as in a noninteracting system. This is because, while levels usually 
have to repel each other except for accidental degeneracies in the non-integrable case, there 
exist enough symmetries to characterize a state that allows the levels to cross in the inte- 
grable case. This should also apply to many-body systems in integrable or nonintegrable 
cases. For instance, Di Stasio and ZotosE have obtained level statistics for spinless fermions 
in one dimension with nearest and next-nearest-neighbor interaction. Thus it is intriguing 
whether the extended spectral flow analysis, which has been shown to be useful in integrable 
systems from the Bethe-Ansatz analysis, could be applied to non-integrable systems as well. 
We find here that the spectral flow is smooth in solvable cases as opposed to wiggly flows in 
nonintegrable cases, but that its period does not change across the integrable-nonintegrable 
boundary. 
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For the question (ii) we conclude here that the extended AB method can indeed detect 
a symptom of changes in characters of quantum states even for non-integrable system. For 
the ID t-J model, the phase diagram has been obtained by Ogata et a/.0 by combining 
the Tomonaga-Luttinger (TL) theory&0 and the Lanczos diagonalization of finite systems. 
With the same method, Ogata et a/.@ obtained the phase diagram of ID t-J-J' model, 
and Sano and Ono0, and independently Kuroki et aS, obtained the phase diagram of 
extended Hubbard model. For the t-J ladder, several authors have obtained the phase 
diagram using the TL theory assuming that the ladder system may be described by a 
ID effective theory&ll The phase diagrams obtained from the extended AB results are 
qualitatively consistent with those predicted by TL theory. 

The organization of the paper is as follows. In section 2 we describe the numerical 
implementation for keeping track of the extended flow. We shall then discuss the result for 
the ID t-J model in section 3, where we perform the Bethe-ansatz analysis to discuss the 
origin of the extended periodicity. Specifically, we discuss the behavior of the extended AB 
period in the non-integrable region from a numerical analysis. Although the period does 
not change across the integrable-nonintegrable boundary, their difference does appear in line 
shapes of the spectral flow. We discuss the results for the t-J ladder (Section 4), ID t-J-J' 
model (Section 5) and extended Hubbard model (Section 6). 



II. NUMERICAL ALGORITHM 

In general we want to treat non-integrable systems. For instance the t-J model is non- 
integrable even in ID except at special points, i.e., J = (U = oo Hubbard model) or 
J = 2t (the supersymmetric t-J model). Thus we must numerically implement the method 
to obtain the extended AB period. Here we recapitulate the numerical algorithm proposed 
in Ref.12, which consists of successive estimations of the energy and wave function. 

1. First solve if ($ = 0)0 = E(<& = 0)0 by a conventional diagonalization technique, 
(e.g., the Lanczos method). 
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2. Predict E n=0 {<5> + A$) by E n=0 = (</>($)\H($ + A$)|0($)>. 

3. Determine ra+ i($ + A<£>) by the power method, 

0($) ^0 n+1 ($ + A$) 

4. Determine E n+ i($ + A$) by + A$)0 n+ i($ + A$) = £ n+ i($ + A$)0 n+ i($ + A$). 

5. If a convergence criterion, |i? n+ i($ + A$) — E n (<$> + A$)| < e is met with typically 
e ~ 10~ 8 , then increase <3> and go to 2. Otherwise, increase n and go to 3. 

Usually the energy converges faster than the wavefunction, so that the above algorithm 
works well. Typically, a few tens of m and n give a convergence. 

Around level crossings a care must be taken. There are two types of level crossings. One 
type occurs between levels that have different symmetries. In this case we can readily go 
across the crossing, since our method looks at the wavefunction, which is a representation 
of all of the symmetries of the system. The second type occurs at a critical point, at which 
a gap opens in the spectral flow. Thus we must question how a crossing turns into an 
anti-crossing from crossing at a critical point. We can systematically keep track of this by 
varying interaction parameters around the critical point. 
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III. ID T-J MODEL 

We first analyze the extended AB period of the ID t-J ring, where the main interest 
is whether a change in the period detects a symptom in characters of quantum states. 
The phase diagram of this model has been investigated by the Tomonaga-Luttinger theory 
combined with numerical results,^ with which we can compare our result. 

The t-J model belongs to a class of the strongly correlated model in which the electron 
repulsion is assumed to be infinite so that no double occupancies of electrons are allowed. 
The Hamiltonian is given by 
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H = -t £ (cl 



c j(T + H.c.) + Jj2(Si ■ Sj 




where c ia creates at the i-th site with spin a and Si is the spin operator defined by Si = 
hcl^a 'aa'Cia' with the Pauli matrix <T aa >. Hereafter we set t — 1 and the phase diagram will 
be parameterized by J and the density of electrons. 

Figure |3], f|, || shows the numerical result for the spectral flow of the ID t-J ring with 
various values of J. Here each system has 4 electrons, while the size of the system is varied 
from 8 (FigfJ), 10 (Fig.[|) to 12 (Fig.^). We can immediately notice three features!!: 

• The extended period, which is originally N a § at J = 0, is abruptly halved as soon as 
J 7^ is turned on. 

• At a certain Jq. > 2, whose value depends on the band filling, a level anticrossing 
appears in the flow line. For 8- and 12-site systems, this makes the period further 
halved into a 1/4 period compared to that for J = 0. In the case of a 10(^ 4,mod4) 
site system, the appearance of the 1/4 period is less perfect. 

• Solvable points (J = 0, 2) are distinguished by smooth, cosine-like flow lines. This 
sharply contrasts with the flow lines in the non-integrable case, where the lines wiggle 
due to ubiquitous level anticrossings. Despite these, the extended AB period does re- 
main the same across the integrable-nonintegrable boundary, which is explicitly shown 
by the fact that the curvature of the anticrossing shows no singularities for all values 
of J considered here. 

In the following sections, we explore why these features should arise. We divide the 
regions into exactly solvable points (J = 0, J = 2), and other regions (normal: < J < 2; 
possibly superconducting or phase-separated: 2 < J). 
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A. J = 



Since the t-J model with J = is equivalent to the U — oo Hubbard model, let us 
start with the ID Hubbard model. Consider an A a -site ring (with N a even for simplicity) 
containing N electrons. The Hamiltonian is given by 

H = -t^2(c\ a c i+la + H.c.) + U^n^nn 

ia i 

where c ia (c| CT ) annihilates (creates) an electron at the z-th site with spin a, n ia = c\ a c ia , U 
the on-site interaction, and t is the transfer energy taken to be the unit of energy hereafter. 

We can analyze the ID Hubbard model in terms of the exact Bethe-ansatz (BA) solution 
due to Lieb and Wu.H If we thread a magnetic flux $ through the ring, the BA equations 
become! 

n e (^pM)j = g e ( 2(A °- A ^) . (2) 

Here M(< A/2) is the number of down-spin electrons, e(x) = (x + i)/(x — i), and kj (X a ) 
is the charge (spin) rapidities, which are roughly a quasimomenta for the charge (spin). 
When U = oo, the charge part of the BA equation (eq.(3.1)) becomes^ 

kjN a = 2izlj + 2tt$/$ + P H , 

M 

P H = ^2tan- 1 (2A a ), 

a=l 

with A Q = 2X a /U. This shows that in the strong-interaction limit the allowed values of k/s 
are quantized in units of 2n/N a . In other words, charge degrees of freedom act like spinless 
fermions. Hence the flow line has a periodicity of A a $ - 

B. J = 2, the supersymmetric point 

J = 2 is special in that the model is integrable with the Bethe-ansatz solution, which 
we first briefly summarize. Laiill and later Schlottmann,il proved that solving the problem 
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at J = 2 reduces to solving a set of (N + M)-coupled algebraic equations for iV charge 
rapidities Vj and M spin rapidities A a , i.e., 

-pew - n r^ + i/2 - 

" ^ - A a + j/2 = _ TT A /3 - A a + j 

^ Uj- - A a - i/2 }} x kp-ka-i 
with Vj = cot(/cj/2)/2. The ground state is given by a set of complex charge rapidities 
(two-strings) and real spin rapidities, where a spin rapidity lies at the center of each two 
string. 

Sutherland!! found an alternative way of solving the problem. There we treat the holes 
and down-spin electrons as dynamical objects in a background of up-spin electrons. In this 
formalism, BA equations become 

exp(ik a N a ) = -[[ J - — 11 

j=i V a ~ Wj - V 2 p=i V a -Vf3 + l 

l _% w j -v p -i/2 

Here v a = tan(/c Q /2) /2 with a — 1, • • • , M\ = Nh + M involves both spin and charge degrees 
of freedom, while Wj is the hole rapidity that represents charge degrees of freedom with 
j = 1,---,M 2 = Nh with Nh being the number of holes. A particle- hole transformation 
connects Sutherland's representation to that of Lai or SchlottmannSi! 

Sutherland's formalism is more appealing for numerical calculations, since (i) the ground 
state involves only real roots, and (ii) we can directly investigate the behavior of the charge 
rapidities. So we follow Sutherland. If we take the logarithm of Sutherland's BA equations 
we have 

Mi M 2 

N a 9(2v a ) = 2nJ a + £ 0(v Q -v )-"£ 9(2(v a - 

P=l 3=1 

Ah 

2nI J = J2d(2(w J -vp)), 
p=i 

with 8(x) = 2tan" 1 (x). J a is an integer (half odd integer) when M\ — M2 + 1 is even (odd), 
while Ij is an integer (half odd integer) if Mi is even (odd). Here we assume I\ < I2 < ■ ■ ■ < 



Im 2 i J\ < J2 < ■ ■ ■ < J M x - For the ground state, we must choose | J a \'s and |ij|'s as close to 
as possible. 

If we thread a magnetic flux $, the left hand side of the charge part of the BA equation 
becomes 

2-kIj = 2-kIj - 27r$/$ . 

In Fig.|6] we show how the rapidities evolve with $ for the ground state. We can see that the 
rapidities change with a characteristic manner with a periodicity of N a § /2. Conspicuously, 
the charge rapidities, Wj's, respond to the magnetic flux in pairs. The pair diverges to ±00 
at some points, around which the pair sandwiches other rapidity v a . Since each pair of 
charge rapidities exchange their places as they cross from —00 to +00, the charge rapidities 
are shuffled like 

{Wx, W 2 , W 3 , U> 4 , • • •} -> {W 2 , W U ttf 4 , W 3 , ■ ■ ■} 

as one extended period ($ = — > $ = N a § /2) is accomplished. The shuffle is similar to 
the case of the small J limit, as we shall see in the next section. We now look more closely 
at the BA equations. 

1. The charge part of the BA equation 

As 2nlj becomes smaller than — Mitt, W\ and W2 vary from —00 to 00. At the same 
time, V\ varies from —00 to 00. In order to fix the range of tan~ x , we redefine the quantum 
numbers Ij so that {h, J 2 , • • • , I M2 } = {-MJ2, -Mi/2 + 1, -Mi/2 + 2, ■ • ■ , -Mi/2+M 2 -l} 
becomes {M x /2-l, Mi/2, -Mi/2 + 1, • • • , -M!/2 + M 2 -2}. When this process is repeated, 
{Ij} comes back to {Ij} when $ = N a § /2. 

2. The spin part of the BA equation 

When v\ varies from —00 to 00, the quantum numbers will be redefined so that 
{Jii J2, ■ ■ ■ , Jah} becomes { J\ + M x — 1, J 2 — 1, J3 — 1, ■ ■ ■ , Jm ± — 1} for the present choice 
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of the Riemann plane. As a set, they do not change. 

Accordingly the BA equations have a periodicity of N a Q Q /2, so that the rapidities return 
to the original values as a set when $ = N a $o/2. Note that in this case, the spectral flow 
corresponds to a trajectory of the holon-antiholon excitations considered by Bares et a/.S 



C. < J < 2 

Remarkably, above halving of the extended AB period for a sizeable J = 2 is seen to 
emerge already for an infinitesimal J from the numerical result. In other words, the t-J 
model with a finite J and the Hubbard model have different extended AB periods. Here we 
show why an infinitesimal J is enough to halve the extended AB period. For a small J/t we 
can treat the J-term perturbatively. One can start from the iV-particle state expanded as 

l i? ) = E-EE-E /<n->i • • II C U • • ■ cUJo) . 

X\ Xjy G\ 

For J = 0, it is known0 that the coefficient / factorizes into charge and spin parts as 

f ai -a N {xi ■ ■■xn) = (-l) Q det[exp(ikiX Qj )]$(yi ■ ■ ■ y M ) , 

where xq\ < xq 2 < xq 3 ■ ■ ■ < xqn are the coordinates of all the electrons with Q being a 
permutation and kj the momentum of a free spinless electron, while y\ < Di - ■ ■ < Dm are the 
coordinates of down spins in the spin configuration. The determinant depends only on the 
positions of particles (xqi, • ■ ■ ,xqn) and not on the positions of down spins (yi, ■ ■ ■ , yM)- 
The spin part &{y\ - • - yM) is identical to the wavefunction of the ID Heisenberg ring of 
length N. 

When a small J is turned on, characteristic change occurs at a special point in the flow, 
i.e., highly degenerate points. Effect of J first appears as a lifting of this degeneracy, which 
can be examined by the first-order degenerate perturbation, as is done@ for the case of the 
ID Hubbard model for U ~ 0. The effective Hamiltonian in the first-order perturbation 
becomesii 
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H c s = -J(n\Y,nin i+1 \m) S fY, S 3 ' S i+i' ( 3 ) 

* i 

where the number operator of spinless fermions, and (n\ ■ ■ ■ \m)sF is the matrix 

element between wave functions of non-interacting spinless fermions having the same energy. 
The summation over j is taken over the compressed spin chain (ignoring the hole sites). If 
we Fourier transform the first part of the effective Hamiltonian as 

i iVa k,k',q 

we can see that when J is turned on, the degenerate states having {■ • • , k, ■ ■ ■ , k', ■ ■ •} can 
mix to give an anti-crossing with the states having {• • • , k — q, ■ ■ ■ , k' + q, ■ • •}. For these 
states to have the same energy, we have (k + k')/2 = ±7r/2(mod27r) = 7r/2 or 37r/2. 

We can compare in Figs.^ ~ ^] the numerical result for J = 0.1 with that for J = 
obtained by the Bethe-Ansatz assuming q = n. It can be seen that the spectral flow at 
J = 0.1 delineates the lower envelope of the J = flows, suggesting that we have indeed 
q = 7r. For example, let us consider the case of 4 electrons in 12 sites (Fig||). When J is 
turned on, the rapidities {kj} change as 

71 37T 57T 771 71 3tT 5n — 7tT 7lT __ — 57T 

*12' 12 ' 12 ' 12 ' *12' 12 '12 ~ ^ ~ ~12~' 12 + ^ = HT' 
at $ = 2.5$ , and 

57T 7lT — 37T — 7T 57T — 7lT 7lT __ — 57T — 37T — 7T 

* 12' 12' ~12~' "12"-' ^12 ~ 71 ~ ~12~' 12 +7T = ~T2~' HP ~Y2* 
at $ = 4.5$ - 

We can indeed give a perturbational argument why the umklapp (q = n) is selected. 

7T 37T 57T TtT 

For J = 0, the state |— , — , — , — ) = |0) is degenerate at $ = 2.5$o with the states 

„ JL — JL — JL — JL — 

7T 37T 57T 77T 37T 4tT 57T 67T — 4lT —571 

— , — , q, h q) = \q), where q is one of vectors, — , — , — , — , or . 

1 12' 12' 12 y ' 12 Hl Hh H '6 6 6 6 6 6 

The matrix for J^niKi+i appearing in eq.(§) spanned by these states are 
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g = 


6tt/6 


57r/6 


4tt/6 


Qtt/Q 


-4vr/6 


-5tt/6 


q = 


( 1 


a~ 6 


a" 5 


a" 4 


a" 3 


a 4 




67r/6 


a 6 


1 


a 


a 2 


a 3 


a 10 


a 11 


5tt/6 


a 5 


a- 1 


1 


a 


a 2 


a 9 


a 10 


4tt/6 


a 4 


a" 2 


a" 1 


1 


a 


a 8 


a 9 


3tt/6 


a 3 


a" 3 


a" 2 


a- 1 


1 


a 7 


a 8 


-4tt/6 


a" 4 


a" 10 


a" 9 


a~ 8 


a" 7 


1 


a 


-5tt/6 




a- 11 


a- 10 


a' 9 


a" 8 


a- 1 


1 / 



where a = exp(i7r/6). This matrix has one nonzero eigenvalue, 7, with the eigenvector 
|m) = (1, a 6 , a 5 , a 4 , a 3 , a -4 , a~ 5 ), while all the other eigenvalues are zero. Hence, when a 
small —JJ2i n i n i+i is turned on, the degeneracy is partially removed with a single level 
peeled off (Fig.|7|). For the flows to be connected continuously, the state |0) must transform 
as |0) — ► \m) — ► \q — it), which exactly implies q = n. The above argument holds for 
any numbers of electrons and sites. Interestingly, this behavior is similar to that for the 
attractive Hubbard model, in which q = tc holds as well. 

The motion or the 'selection rule' for the rapidities can be seen more directly in Fig|| 
where two fe's pair off and respond to $ hand in hand. We can interpret this as an effect 
of J working as an attractive interaction. As a consequence, a finite J halves the period of 
the spectral flow. Namely, the fc's for the ground state change from {k±, k 2 , k 3 , fc 4 , ■ • • , } to 
{k 2 , h, h, k 3 , ■ ■ ■ , } for $ = -> N a %/2. 

This behavior is reminiscent of the charge rapidities in Fig.6, which can be defined when 
the system is integrable at J = 2, the supersymmetric case. For J = 2, Bares et a/.El have 
termed the phase as 'weakly bound electron pairs' having antiparallel spins because charge 
rapidities by Lai's representation form two strings. In terms of the extended AB period, the 
numerical result shows that the period is the same over the whole region of < J < Jc- In 
this sense the 'weakly bound pairs' region extends over < J < Jc- 
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D. J > 2 



The present result, Figs.| ~ exhibits a period halving at a critical Jci> 2). In order to 
determine Jc, we have evaluated the energy gap A in the following way. First, we identify, 
with a linear extrapolation, the point Eq at which levels would cross if the level repulsion 
were absent. The energy gap is then estimated as A = E — Ex, where E\ is the energy of 
the repelled level. We can next plot A as a function of J and determine Jc as A( Jc) = by 
extrapolation. The results are Jc — 2.30 for 4 electrons in 8 sites, Jc — 2.22 for 4 electrons 
in 10 sites, Jc — 2.14 for 4 electrons in 12 sites. 

According to the phase diagram obtained by Ogata et the superconducting correla- 
tion function becomes dominant (K p > 1) in the region J > J c , where J c > 2 is a function 
of the band filling. The phase boundary is identified as the trajectory of K p = 1 in terms of 
the exponent K p that characterizes the TL liquid. The boundary between superconducting 
and phase separated regions is more or less pararell to the trajectory of K p = 1. 

Figure |9] superposes the period halving points on the phase diagram (trajectory of K p = 
1) the TL resultEl. We can see the halving points fall on a line near the trajectory of K p = 1, 
which is pararell to the boundary line between superconducting and phase separated regions. 



IV. T-J LADDER MODEL 

We now turn to the t-J ladder model (for a review of the ladder, see, e.g., Ref.37). The 
original motivation to introduce ladders came from an expectation that a gap in spin ex- 
citations may be formed in such systems, leading possibly superconductivity when carriers 
are doped. The importance of the spin gap has in turn conceived for the high-T c cuprates. 
Specifically, the t-J ladder has been predicted to have superconductivity associated with a 
spin gap when there are even number of chainsHH The interests are heightened when a 
class of cuprates are found to possess the ladder structure in the Cu-0 network.0^ii Very re- 
cently, a superconductivity has been detected in (Ca,Sr)i 4 Cu2404i under high pressures.3^1 
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This has kicked off numerical works with the exact diagonalization or the density- 
matrix renormalization group (DMRG) method-lllMBEl^El Since the system is quasi-lD, 
the Tomonaga-Luttinger (TL) liquid theory has also been applied.@~0 Some results have 
predicted existence of a superconducting phase, in which a spin gap is also shown to survive 
the doping near the half filling.0 However, the ladder is not purely ID, the validity of the 
TL approach has to be justified. Hence, another method detecting the pairing is desirable. 

Thus it is most intriguing how the extended AB period test predict the transition for 
the t-J ladder. Kusakabe and Aoki@ obtained a preliminary result that the extended AB 
period changes around the phase boundary predicted by Sano@ We discuss the extended 
AB period of this model in more detail. 

The t-J ladder is defined, in standard notations, by the Hamiltonian, 

H=~tJ2 P G (ct ai(7 C i+ i iQ)CT + h.C.)P G ~t±J2 P G(4,l,a C i,2,c + h.C.)P G 
i,a,a i,tr 

+ J ^^(S-t,Q ' Sj_|_l jQ ~1T-i+l,a n i,a) 
i,a 

+ J± XX S '.i " S *.2 - ^n iA n ii2 ), (4) 

i 

where a labels two legs of the ladder, while i labeling the rung runs from 1 to N a . Doubly 
occupied sites are totally excluded by the Gutzwiller projection Pq as in the single t-J 
model. 

Here we take the same intra- and inter-chain electron transfer, t = t±, and the same intra- 
and inter-chain superexchange interaction, J — J±, for simplicity. The AB flux threaded to 
the ladder wound along the length introduces a Peierls phase as 

t (> \ t 

C i,a,a C i+l,a,a ex P ^jy ^ J C i,a,cr C i+l,a,a ■ 

We have numerically obtained the extended spectral flow for finite t-J ladders. Results 
depends on the number of electrons. 
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A. Two-Electron Systems 



Result for two electrons in an N a x 2 = 8 x 2-site system is shown in Fig.|l(] We look 
at the J-dependence of the flow of the ground state, starting from the t model (the system 
with no superexchange interactions, J = 0, but the double occupancies still excluded). For 
J = the flow line has a periodicity of N a $ , which is the period expected for individually 
moving particles. Remarkably, a small J(= OAt) is enough to halve the periodicity. This 
behavior, which is common to two-electron systems on In x 2-site ladder in: an integer), 
is caused again by a pair of level anticrossings in the flow. The change in the period from 
iV a $o to iV a $o/2 is similar to the result for the single-chain t-J model (Section 3), in which 
the period is also halved for an infinitesimal J. For general numbers of electrons in a ladder, 
however, the period N a Q is not always observed even for J = 0. For some numbers of 
electrons, the halved period, N a Q /2, is often detected (see below) for the t ladder, although 
t model has no explicit source of attractions between electrons. 

B. Four-Electron Systems 

We move on to the spectral flows for four electrons for an 8 x 2-site system (1/8-filling) 
in Fig.|TT] for various values of J. We can see that there exist a pair of level anti-crossings 
at Jc(— 2.6t), which reduces the extended period into iV a $ /4. 

For a 6 x 2-site system (1/6-filling) in Fig.|l2|, a similar appearance of a pair of level 
anticrossings occurs at Jc — 2.04t,@ although the period halving is less ideal. This can be 
regarded as a pseudo 1/4-periodicity, in which the exact 1/4 period is inhibited due to the 
incommensurability of iV a with 4. 

In Fig.[13] we compare the Jc obtained from the extended period halving with the phase 
diagram obtained from the TL analysis by Hayward et a/.H and independently by Sanc@. 
The value of Jc falls on a line near the trajectory of K p = 1 and nearly pararell to the 
boundary between superconducting and phase-separation regions. 
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Here we should comment the following. The period is already halved at J = (the t 
ladder) at least for the 6x2 system.S This should be related to the Hubbard ladder (the 
Hubbard model on a double chain), since the t model is identical to the Hubbard model 
with U — oo. Thus the Hubbard ladder should have the halved period as well for U — oo, 
while for {7 = the extended period is iV a $ - This implies that a critical point has to exist 
at which the extended period changes in the Hubbard ladder. A preliminary calculation 
for a 6 x 2-site Hubbard ladder with 4 electrons indeed indicates a halving of the extended 
period at a finite U c . If this is true, we have another example in which the property of the 
ladder is distinct from the single chain, since the period is always iV a $o i n the single-chain 
Hubbard model with U > 0. 



C. Quarter filling 

At the quarter filling (N a electrons on a N a x 2 ladder), a peculiar behavior appears, 
i.e., the period becomes the smallest possible $o irrespective of J and the systems size 
(4 x 2, 6 x 2 and 8x2 sites). As an example we show the ground-state spectral flow for 
a 6 x 2-site system with 6 electrons in Fig.|l4]. The period $o is what is expected for the 
charge-gapped system (e.g., a Mott insulator). In this sense the $o period may come from 
the commensulability of the filling, although we could not find any characteristic structure 
in the CDW correlation (n^n^) , which is almost flat in these small systems. The period 
remains the same when we modify the interaction parameters in the region t < t±, J < J±. 
On the other hand, if we approach to the independent chains in the region t > t±, J > J±, 
a branch having a longer period comes downwards to cross the lowest $o _ P er i°d line. This 
is natural because in the limit of t ^> t±, J 3> J±, the extended period for the ground state 
should be that for the ID ring with N/2 electrons. 
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V. T-J-J' MODEL 

Ogata et al have studied ID t-J-J' model,!! where they added a frustrating next nearest- 
neighbor exchange interaction J' to the t-J model. Like ladder systems, this has been 
motivated from an expectation that we can tune the spin gap by introducing J'. J' is in fact 
shown to give rise to a spin gap at half-filling, which survives for small doping and arbitrary 
values of J. In this section, we apply our method to t-J-J' model to discuss whether the 
extended AB period can detect a symptom of changes in characters of quantum state. 

In particular, the halving point, Jc, of the extended AB period from N a § /2 to N a ^ Q /A 
is of interest. (As in ID t-J model, the extended AB period halves to N a &o/2 right after the 
J and J' set in.) We fix J'/ J = 0.5, because Ogata et a@ determined the phase diagram 
for this case. The present result shows that Jc = 3.4t for 4 electrons in 8 sites, Jc = 2.8t for 



4 electrons in 10 sites, and Jc = 2.6i for 4 electrons in 12 sites. In Fig.[T5|, we superpose the 
present result on the phase diagram due to Ogata et a@. We can see that Jc determined 
by the extended AB effect falls near the neighborhood of the boundary of phase separation 
obtained by the exact diagonalization method. 



VI. EXTENDED HUBBARD MODEL 

So far we have studied models for which the phase diagram contains spin-gapped regions. 
However, we have not studied the behavior of the AB period in such regions so far. It is 
interesting to investigate how the extended AB period behaves in the region where spin gap 
opens, since the opening of a spin gap is regarded as a tendency toward attraction between 
electrons. In the extended Hubbard model, where we consider the off-site (nearest-neighbor) 
interaction V in addition to the on-site Hubbard U, the phase diagram has both spin-gapped 
and spin-gapless regions bisected by the so-called Luther-Emery line. Thus the model is ideal 
for studying how the existence or otherwise of the spin gap affects the extended AB effect. 

The extended Hubbard model is defined by the Hamiltonian, 
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H = -t ( c L c jv + H.c.) + U 53 n il n H + V 53 n «^+i<x- 

<i,j>cr i icr 

According to the weak-coupling theory, the superconducting correlation dominates in the 
region 

U + 2[2- cos{2k F ))V < 
and the spin gap opens in the region 

U + 2V cos(2A; F ) < 

For the strong-coupling region, the phase diagram has to be obtained numerically. Mak- 
ing use of the method same as the ID t-J study by Ogata et a£l, Sano and Onoll, and 
independently three of the present authors!!, have numerically determined the phase dia- 
gram. 

Here we focus on the region with V < for which the pairing is expected. When V set in 
the extended AB period halves to N a Q /2 regardless of the value of U. We have determined 
the point where the extended AB period is halved from N a Q /2 to NJ&q/A for 4 electrons in 
8 sites, 4 electrons in 10 sites, and 4 electrons in 12 sites. In Fig.[L6|, we superpose our results 
on the phase diagram due to Sano and Onoil. The halving point of the extended AB period 
is seen to run roughly parallel to the normal-superconductivity and superconductivity-phase 
separation boundaries, and falls on the superconductivity region.0 

VII. CONCLUDING REMARKS 

In four types of strongly correlated electron systems in ID, namely 

1. t-J model, 

2. t-J ladder, 

3. t-J- J' model, 
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4. extended Hubbard model, 

we have obtained the extended AB period both from the Bethe-ansatz analysis where pos- 
sible, and from numerical calculations. While the extended AB period is the system length 
times the flux quantum for noninteracting systems, we have found the existence of the 
boundary across which the period is halved or another boundary into an even shorter period 
on the phase diagram for these models. If we compare this result with the phase diagram 
predicted from the Tomonaga-Luttinger theory, the halved period (or shorter periods) has 
a one-to-one correspondence to the pairing (phase separation or metal-insulator transition) 
in these models. Namely, when the halved (quarter) AB period is detected, the correspond- 
ing TL phase diagram has a paired (phase-separated) region. Such an applicability of the 
extended AB test extends to the cases where the conventional anomalous flux quantization 
gives ambiguous results as is the case with the Hubbard or t — J models. 
Surprising aspects of these results are: 

1. The extended AB flows involve high-energy states, while the TL theory is only intended 
for the low-energy physics. 

2. The extended AB period does not change across the integrable-nonintegrable boundary 
despite the totally different level statistics. 

An important future problem is how the extended AB analysis can be coupled with 
Berry 's geometrical phase. Korepin and have in fact calculated for the XXZ Heisenbere; 
magnet Berry's phase from the Bethe-ansatz solution based on the periodicity ($ = 2$o) 
found by Sutherland and Shastry. Berry's phase, 7, is expressed in a standard way as 



We can in principle calculate Berry's phase for electron systems such as the Hubbard model 
similarly. The geometrical phase will help identifying the nature of the many-body states in 
a gauge field theoretic manner. 




(6) 
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FIGURES 



FIG. 1. An Aharonov-Bohm flux pierced in an opening of a ring. 

FIG. 2. The extended AB flow for a 10-site Hubbard ring with 6 electrons 12 . The flow lines 
represent a repulsive interaction (U > 0), an attractive interaction (U < 0), and the noninteracting 
case. 

FIG. 3. Numerical result for the spectral flow for a ID t-J ring with 4 electrons in 8 sites. 
FIG. 4. Numerical result for the spectral flow for a ID t-J ring with 4 electrons in 10 sites. 

FIG. 5. Numerical result for the spectral flow for a ID t-J ring with 4 electrons in 12 sites. 

FIG. 6. The rapidities for the hole charge {wj\ thin lines) and for spin-charge degrees of freedom 
(v a : bold lines) against the flux for the ground state for a ID t-J ring with 4 electrons in 10 sites. 
Tiny gaps in the non-diverging branches correspond to the points where a pair of branches diverge 
to ±oo. 

FIG. 7. The change of the spectral flow when small J is turned on. 

FIG. 8. The motion of charge rapidities for a finite J in ID t-J model. 

FIG. 9. The present result for the halving points (open circles) of the extended AB period in 
ID t-J model is superposed on the phase diagram obtained by Ogata et al. 18 

FIG. 10. The extended AB flow for an 8 x 2-site t-J ladder with two electrons for J = 0, 
J = OAt and J = 0.8t, respectively. 

FIG. 11. The extended AB flow for an 8 x 2-site t-J ladder with four electrons for 
J jt = 0, 2.0, 2.6, 3.0, 3.4 from top to bottom. 
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FIG. 12. The extended AB flow for an 6 x 2-site t-J ladder with four electrons for 
J = 0, 2.0, 2.2, 2.6, • • • , 3.8 from top to bottom 

FIG. 13. Halving points, Jc, of the extended AB period (big crosses) in the t-J ladder are 
compared with the trajectory of K p = 1 (the normal-superconducting boundary) obtained by the 
exact diagonalization (solid squares by Hayward and Poilblanc 25 , solid triangles by Sano 26 ). Phase 
separation boundary (K p = oo) are also displayed by open squares or open triangles. 

FIG. 14. The extended AB flow for an 6 x 2-sites t-J ladder with six electrons (the quar- 
ter-filling) for J = and J = t. The period is <J>o irrespective of J. 

FIG. 15. The present result for the halving points (filled circles) of the extended AB period in 
ID t-J- J' model is superposed on the phase diagram obtained by Ogata et al. 22 

FIG. 16. The halving point of the extended AB period in the extended Hubbard model for 
four electrons in an 8-site system (band filling n = 0.5; open circles), 10-site system (n = 0.4; 
open triangle), and 12-site system (n = 0.333; open square) is superposed on the phase diagram 
obtained by Sano et al. 23 and Kuroki et al. 24 for n = 0.5. For n = 0.5, The line U = is the 
boundary between the spin-gapped and spin-gapless regions. 
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